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Abstract 

The generating functionals for the local composite operators, $ 2 (x) 
and $ 4 (x), are used to study excitations in the scalar quantum field 
theory with \Q 4 interaction. The effective action for the composite 
operators is obtained as a series in the Planck constant h, and the two- 
and four-particle propagators are derived. The numerical results are 
studied in the space-time of one dimension, when the theory is equiv- 
alent to the quantum mechanics of an anharmonic oscillator. The 
effective potential and the poles of the composite propagators are ob- 
tained as series in h, with effective mass and coupling determined by 
non-perturbative gap equations. This provides a systematic approx- 
imation method for the ground state energy, and for the second and 
fourth excitations. The results show quick convergence to the exact 
values, better than that obtained without including the operator <5 . 



1 Introduction 

The quantum mechanical anharmonic oscillator (AO) can be regarded as a 
quantum field theory of a real scalar field with a classical action given by 



when the dimension of the Euclidean space-time n is equal to one. For 
this reason the AO is frequently used as a testing ground for approximation 
methods in quantum field theory. The field-theoretical perturbation method 
enables us to calculate Green's functions and their generating functionals 
to an arbitrary order in the coupling constant A. Within this approach the 
perturbative expansion of the ground state energy of the AO is easily obtained 
by calculating the effective potential in the loop expansion, but a derivation 
of excitation energies from zero modes of the appropriate Green's functions 
is not so straightforward Q, in contrast to the Schrodinger approach, where 
perturbative calculations for excited levels are of the same difficulty as for 
the ground state. 

For a field-theoretical investigation of excited modes, the N— particle 
functionals which depend on some artificial sources coupled to composite 
operators, $(xi)$(x 2 )...$(x7v), have been introduced [0, |]. These function- 
als generate composite propagators through differentiation with respect to 
the sources. The simplest Green's functions, with the simplest relations to 
multi-particle eigenmodes, are generated by the effective action for composite 
operators, defined through a Legendre transform of the connected generat- 
ing functional. For N < 4 the effective action has been shown [^, |]] to be 
a sum of A-particle-irreducible vacuum diagrams with the exact propaga- 
tor. In the case N = 2, studied in the classic paper of Cornwall, Jackiw 
and Tomboulis ||, the Gaussian effective action (and the Gaussian effec- 
tive potential ||) is obtained as an approximation of the two-loop result. It 
is, however, difficult to discuss a full content of the two- loop result, and to 
proceed to higher operators within such an approach because of the highly 
non-trivial integral gap equations for vacuum expectation values of composite 
fields. 

The local N— particle operators $ (x) have been also considered 0; how- 
ever, diagrammatic rules for the effective action were lacking, since the Leg- 
endre transform cannot be performed explicitly in this case. Only recently, 
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an expansion of the effective action for the operator $ 2 (x) in terms of two- 
particle-point-irreducible (2PPI) diagrams was given || |9|. The result is 
implicit but enables us to calculate the effective potential || and the two- 
particle composite propagator |IJ. The Gaussian effective action is obtained 
already in the one-loop approximation, and the calculation of post-Gaussian 
corrections is easier than for the bilocal operator, since the gap equation for 
the vacuum expectation value of the local composite field is algebraic. Un- 
fortunately this approach cannot be extended to N > 2, and all one can do 
in that case is to perform the inversion in the Legendre transform order by 
order in the Planck constant ft. For multilocal operators this uphill task can 
be done explicitly, and the well-known expansions in N— particle-irreducible 
diagramms are recovered [IT]. For the local composite operator $ 2 (x) the 



method appears implicit; nevertheless it makes a discussion of many applica- 
tions possible [fr2] , jT3[ . Even the diagrammatic rules for the effective action 
(different from that of the 2PPI expansion) have been established [II], [14] . 
The aim of our paper is to extend the inversion method to include the com- 
posite operator $ 4 (a;). 

In Section [|, after a brief exposition of the generating functionals for 
composite operators, a diagrammatic representation of the effective action 
for the <& 2 {x) operator is obtained explicitly to the order ft 5 , and the (inverse) 
two-particle composite propagator is derived up to ft 4 . Then we calculate the 
effective action for two operators, $ 2 (a;) and $ 4 (x), to the order ft 4 . This 
enables us to obtain the two- and four-particle propagators up to ft 3 and ft 2 , 
respectively. Having in mind an application of the formalism in relativistic 
quantum field theory, we discuss the scalar theory in the space-time of n- 
dimensions, keeping n arbitrary as long as possible. It is only in Section |3], 
when calculating excitation energies, that we set n — 1 to consider the AO. 
The poles of the composite propagators are determined in successive orders 
of ft and the resulting energies are compared with the exact spectrum of 
the oscillator. The approximations for the ground state energy and for the 
second excitation, obtained from the effective action for $ 2 (x), are in good 
agreement with the exact results; however, in this approach there is no simple 
way to derive the fourth excitation. The advantage of using the effective 
action for two operators, $ 2 (x) and $ 4 (x), lies in the fact that it generates 
directly the four-particle propagator, providing an excellent approximation to 
the fourth excitation energy of the AO. Moreover, the approximations to the 
ground state and to the second excitation obtained from the effective action 
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for the operators $ 2 (x) and $ 4 (x) are considerably better than that obtained 
from the effective action for the operator $ 2 (x) only. Our conclusions are 
summarized in Section [|. 



2 The effective action for local composite op- 
erators through the inversion method 

2.1 Generating functionals for composite operators 

The vacuum functional for the local composite operators $ 2 (x) and <& 4 (x) 
can be represented by a path integral 

Z[J,K] = e ^ w[J ' K] = J D <& e l[-m^p{x)*\*)d n x+±jK{ x )z\ x )dr x ]_ ^ 

The connected generating functional, W[J, K], is given by a sum of vac- 
uum diagrams in configuration space, obtained with an inverse propagator 
Gj (x,y) = (— <9 2 + m 2 — J(x))8(x — y) and a four-point vertex given by 
{K(x) — 24X)5(x — y)8(x — w)8(x — z); in Fig. ?? we show the diagrammatic 
expansion to the order h 5 . For simplicity, we consider here the case of unbro- 
ken reflection symmetry, when the vacuum expectation value of the scalar 
field vanishes; this is why we do not introduce a current coupled to an ele- 
mentary field <3?(x). This, however, will limit our investigation to the even 
excitations of the system. 

The effective action is defined through a Legendre transform 



r[A, A] = W[J,K] - - Jj(x)A(x)d n x- — J K(x)k(x)d n x 

*K{x)A 2 (x)d n x, (3) 



where the background fields are defined as 
5W 

hA(x) = 2—— =< $ 2 (x) >j :K 
dJ{x) 

h 3 A(x) = 2A^^-3h 2 A 2 (x)=<^(x)> JiK -3<^ 2 (x)>l K (4) 
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and < ... >j,k denotes the expectation value in the presence of external 
currents J(x) and K(x). 

The external sources J and K were introduced artificially in order to 
define the background fields (f|) and the effective action as a functional of 
those (||). The physical quantities have to be calculated at J = K = 0, or 
equivalently at the values of the background fields, A and A, for which the 
gap equations 



5T 



5A{x) 



h 

J{x) - —A{x)K{x) = 0, 



5T 

5A(x) 
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(5) 



(6) 



are satisfied. The effective action for the composite operators, T[A,A], de- 
termines the vacuum energy density 



E n 



r[A, A] 



Jd n x 



(7) 



A=const,A=const 



If another solution to the gap equations (A(x) = A + DA(x), A(x) = A + 
DA(x)) exists in the vicinity of (A, A), then by (|5]) and @ the stability 
condition J7|, is obtained in the form 



I 



T 44 (x,y) 



DA(y) 
DA(y), 



dy 



(8) 



The matrix T, given by 
^22(x,y) 



5 2 r 



SA(x)SA(y) 

^42(x,y) 
5 2 T 



n8(j{x) + \K{x)A{x 



A, A 



5 2 T 



5A(x)6A(y) 



5A(y) 
h 3 5K(x) 



A, A 



A, A 



24 8A(y) 



A,A 



5A(x)5A(y) 



h 3 SK(x) 



A, A 



24 5A(y) 



(9) 



A, A 



is an inverse of the composite propagator matrix, defined as 



W(x,y) 



W 22 (x,y) W M (x,y) 
W A2 (x,y) W u (x,y) 



8 2 W 



8J(x)SJ(y) 
S 2 W 



6K(x)SJ(y) 



S 2 W 



SJ(x)SK(y) 
5 2 W 



0.0 

00 SK(x)SK(y) 



0.0 
0,0 
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< T® 2 {x)<S> 2 {y) > con < T$ 2 (x)$ 4 (y) > L 

< T®\x)&{y) > con < T<$>\x)<$>\y) > c 



In Fourier space the stability condition M) is written as 



r 2 2(p) r 24 (p) 



(10) 



£>A(p) 







£>A(p), 








(11) 



hence the zero modes of T{p) are related to the excitations of the system. 
The matrix T(p) can be diagonalised, providing the inverse of the two-particle 
propagator in the form 



and that of the four-particle propagator equal to 

h 3 r 8 (p) = r u ( P ) - r 42 (p)rj 2 1 (p)r 2 4(p). 



(12) 



(13) 



The effective action for the composite operators, T[A, A], provides a con- 
venient tool to study the ground state and excitations. Unfortunately, a sys- 
tematic approximation scheme for the effective action is intricate, although 
the connected generating functional, W[J,K], can be easily obtained in terms 
of Feynman diagrams in the loop expansion. The main difficulty lies in per- 
forming the Legendre transform (|3|), i.e. in eliminating the sources J and K 
in favour of A and A. Here we shall do this order by order in h, representing 
the vacuum expectation values A and A by the series 



A 



fc=0 



(14) 



A 



k=0 



(15) 



where the coefficients A^[J, K] and A^[J, K] can be easily found from Eq. |1 
The diagrammatic representation of A and A to the order h 4 and ft 2 respec- 
tively is given in Fig. ??. The inverted series are written as 



J=Eft*J ( *)[A,A], 

k=0 



(16) 
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K= y £h k K (k) [A,A], (17) 

k=0 

and the coefficients J^[A,A] and K^)[A,A] have to be determined order by 
order in ft. It has to be noted that only the lowest-order inverse functionals, 
Jo [A, A] and K [A,A], are required; higher order coefficients are obtained 
as functionals of them. This establishes an algorithm for calculating the 
effective action to an arbitrary order in h even if Jo [A, A] and Kq[A, A] cannot 
be obtained explicitly, which is the case for local composite operators. 



2.2 The effective action for $ 2 (x), T[A] 

The effective action for the two-particle operator $ 2 (x), T[A], has been re- 



cently examined with great care by Okumura |LT] and Yokojima [[L4[]. As 
an introduction to the case of four-particle operator, discussed in the next 
section, we review here the calculation of T[A] by performing the Legendre 
transform of the connected generating functional W[J] to the order h 5 . In 
the definition of the effective action for the two-particle operator, T[A], only 
the current J{x) coupled to $ 2 (x) is required; therefore K(x) has to be set to 
zero in the formulas for the generating functionals (0-??) and only Eq. (|14]) 
has to be inverted, to determine a functional J [A] to the given order in K. 
Using the identity 

oo oo 

A = £^A (fc) [Y>y w [A]], (18) 

fc=0 i=0 

one obtains an infinite sequence of equations 

A = A (0) (19) 



A' (0) J(i) + A (1) =0, (20) 

l (0)"(2) T -^ (0 )(J(1)) 2 + A (1) , 



A' (0) J (2) + ^A" (J {1) ) 2 + A' (1) J (1) + A {2) = 0, (21) 



A (0 )J(3) + A (0) J ( i)J(2) 

+ ^;V J (i)) 3 + A (i) j(2) + ^ A (i)( J (i)) 2 + A (3) = 0, (22) 
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A'(0)^ (4) + ^(0)(2 t /(i) t / (3) + (^)) 2 ) 

+ \^){J{i)?J{2) + ^ A '(o)(^(i)) 4 + A'(i/ 3) + ^'(dW 

1 

— L 

2 



A (l) J W J m + A (2)^(2) + ^A (2) (J (1) ) 2 + A (3)^(l) + A (4) = (23) 



where A(fc)[Jo] are given in Fig. ??. Here and in the following the space- 
time indices and integration over them are suppressed for notational simplic- 



ity. The coefficients J(j), calculated from Eqs. |20|- £3], can be represented by 
Feynman diagrams in configuration space with the (inverse) free propagator 

Gj^x,y) = (-d + n 2 (x))5(x-y), (24) 

where 

ft 2 [A] = m 2 - J [A]. (25) 

The result to the order h 4 is shown in Fig. ??; the inverse two-particle propa- 
gator, attached to an external point, is that of two free particles of an effective 
mass tt(x). In Fig. ?? we show the effective action T[A] to the order h 5 , ob- 
tained by eliminating J with the use of Fig. ??. In agreement with the rule 
proved by Okumura [JTTJ] , T[A] appears to be a sum of one- vertex- irreducible 
diagrams with the inverse composite propagator attached to an external 
point, and the two-point pseudo- vertices of order i (i > 2) inserted in all 
possible ways. By way of digression, it can be observed that all pseudo- vertex 
insertions in T[A] can be summed by taking a dressed inverse propagator in 
the form G~^ m (x, y) = (-d 2 +m 2 - J (x)- J 2 {x)- J^{x)- J A (x) + ...)5(x-y) = 
(— <9 2 + m 2 — J(x) + Ji(x))5(x — y), and the result coincides with that ob- 
tained in the 2PPI expansion |9| which establishes a relation between the 
two methods. 

The definition of the effective action T[A] is implicit, since the functional 

A[J ](x) = G Jo (x,x), (26) 

obtained from (|19D , cannot be inverted for a space-time dependent Jo- How- 
ever, using the fact that all dependence on J is through J , Eq.ES enables us 



to calculate tx44 

5A(y) 



, and hence the two-particle inverse propagator 

A(x)=A 



T\x-y) = T 22 {x-y) = --fd n z 



h f 5J(x) 5J (z) 



5J (z) 5A(y) 



• (27) 

A(x)=A 
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Since a constant A corresponds to constant Q, the vacuum energy density (0), 
as well as the Fourier transform of the two-particle two-point Green's func- 
tion, T 4 (p), can be represented by the Feynman diagrams in momentum space 
with the propagator 2 _^q2 • The vacuum energy density is given by the same 
set of diagrams as the effective action (Fig. ??), but in momentum space. 
The result for T 4 (p) calculated from fl27|) to the order h 3 is shown in Fig. ??. 
In this case ( |T9"D can be written as 

A = /(l)v + fi 2 ' (28) 

and can be used to eliminate A in favor of Q in all the quantities considered. 
The gap equation (|5]), which becomes algebraic, determines the value of fl 
which corresponds to vanishing external current. 



2.3 The effective action for Q 2 (x) and $ 4 (x) operators, 

r[A,A] 

In order to calculate the effective action for the operators <E> 2 (x) and $ 4 (x), 
r[A,A], both the series given by (|l4l) and (|15|) have to be inverted. Using 
the identities 

oo oo oo 

A = £ n fc A (fc) E tiJ®[A, A], ]T tiK^A, A]]. (29) 

k=0 i=0 k=0 

and 

oo oo oo 

A = ]T ft fe A (fc) E fc' J (0 [A, A] , 2 &%■) [A, A]] . (30) 

k=0 i=0 i=Q 

we obtain two infinite sequences of equations 

A = A {0) (31) 

M(0) -J 1 + A (1) =0, (32) 



6 J, 



o 



TjT J(2) + 2^JT ( J(1)) + (1) + W W + (2) ~ ' (33) 
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and 

A = A ( o) (34) 
^A(o) T SA {0) 

5A (0) _ , 15 2 A (0) 2 5A (0) 5 2 A (0) 5A (1 



(36) 

• ■ • j 

where we made use of the fact that A( ) and -j^j- do not depend on Kq. Tak- 
ing A(fc)[Jo, i^o] and A^[J , K ] from Fig. ??, we calculated the coefficients 
J(i) and from the above equations. The results for J and X to the order 
h 3 and ft 2 respectively are represented in Fig. ?? by Feynman diagrams in 
configuration space, where the composite propagators are those of two or four 
free particles. The effective mass in a free propagator, fl[A] (0), is the same 
as in the case of T[A] considered in Section [2.2| , and the four-point vertex 
is (Kq(x) — 24\)5(x — y)5(x — w)S(x — z). In Fig. ?? we show the effective 
action, T[A, A] to the order ft 4 , which is obtained as a Legendre transform of 
the connected generating functional by using Fig. ?? to eliminate J and K 
in favour of A and A. It seems that, as in the case of T[A], graphical rules 
could be established, but we do not disscuss this point here. 

The expression for T[A, A] given in Fig. ?? is implicit, with the relations 
for J [A,A] and i^ [A,A] given by ( |3l]) and ([34]) respectively. The vacuum 
energy density (0) is represented by the same set of Feynman diagrams as 
the effective action (Fig. ??), but in momentum space with the propagator 
p i+ n i ■ Using the fact that all dependence on sources is through J Q and K , 
the implicit expression for T[A,A] enables us to derive the inverse of the 
composite propagators matrix, defined for constant A and A, and its Fourier 
transform, T{p) can be diagonalized. In Figs. ?? and ?? we show the Feynman 
diagram representations for the inverse of the two-particle propagator, T A (p) 
(up to h 3 ), and for that of the four-particle propagator, T 8 (p) (up to h 2 ), 
respectively. 

For constant A and A, Jo, Kq and Q are also space-time independent, 
thus ( 131]) becomes (28), and ([34]) can be written as 



A = A ( o)(Jo,Ko) 
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(K - 24A) J ■ 



d n pd n qd n r 



( 27r )3n( p 2 + fi 2)( g 2 + fi 2)( r 2 + Q2)((p + g + r )2 + fl 2) 

and inverted to obtain 

A 



(37) 



K = 24A + 



I 



d n pd n qd n r 



(38) 



(27r) 3 "(p 2 +n 2 )(q 2 +^ 2 )(r 2 +^ 2 )((p+g+r) 2 +n 2 ) 



The relations (2"8"|) and (58"|) can be used to eliminate A and i n favor of f2 
and A in the calculated quantities. The values Q and A have to be determined 
from the gap equations (|5]) and (^|), which become algebraic. 



3 The quantum-mechanical anharmonic oscillator 

In the space-time of one dimension the A$ 4 theory is equivalent to the quan- 
tum mechanical anharmonic oscillator with a Hamiltonian given by 

1 1 

H = -p 2 + -mV + Ax 4 . (39) 

The exact element of the composite propagator matrix in the Euclidean for- 
mulation can be represented as 

Wijty-f) = <0|Tx J (t)x J '(t')|0> - <0|x 4 (t)|0><0|^'(t / )|0> 

oo 

= £ <0|x i |/cXA;|^|0>e- |t " t ' |£fe , (40) 

k=l 

where e& = — E$ is the excitation energy of the k— th state, \k 3>. The 
Fourier transform of the above propagator is given by 

W M = E J 2 t p2 — > ( 41 ) 

fc=i p £ fc 

and has an infinite number of poles at imaginary momenta whose absolute 
values determine all excitations of the system. Equivalently the full spectrum 
can be obtained by looking for zero modes of the exact inverse propagator 
matrix, T{p). However, if we are working with an approximate propagator, 
which may have only a finite number of poles, this provides an approximation 
to some part of the energy spectrum only. Different, but usually better, 
approximations can be obtained from zero modes of the T matrix, and such 
an approach will be used in our work. 
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3.1 The approximations of the AO spectrum 
3.1.1 The approximations derived from T[A] 

The vacuum energy density of the scalar field in the one dimensional space- 
time gives the ground state energy of the AO. Using the effective action for 
the operator $ 2 , T[A], given in Fig. ?? the ground state energy of the AO is 
calculated to be 

„ n ,„ 2 2 ^A „,«. 2a2 3A 2 ft 3 27A 3 ft 4 1923A 4 ft 5 

E = h (tt 2 -m 2 )h— + 3Xh 2 A 2 - + + - (42) 

2 V J 2 80 5 160 8 1280 11 V ; 

where A is space-time independent, and by Eq. (p8| is given by 

A - h < 43 > 

The two-particle inverse propagator, from Fig. ??, becomes equal to 

n fi / 2 l2t2 160O 2 + p 2 „ „ 3 71680 4 + 1760V + p 4 

" r < P » = 2 < P + 4 " » + 6Afi " 3A R V(16?^) + 27A ?l 2n.(lOT + £)» 
4 4 5198561280O 8 + 4167464960V + 118792320V + 1606920 V + 64 lp 8 

160 9 (160 2 +p 2 ) 3 (360 2 +p 2 ) + """ ( ^ 

and its zero determined to the order fr 4 results in an approximation to the 
second excitation given by 



_ 3hX 87ft 2 A 2 2547ft 3 A 3 40169m 4 A 4 
€2 = 2fi+ 0^-^ + ^ 640^ + - (45) 



The effective mass O has to be determined from the gap equation ([5]) which, 
after using fl43|) , in the one-dimensional space-time becomes 

320 12 - 32m 2 10 - 192Aft0 9 + 240A 2 ft 2 O 6 - 1728A 3 ft 3 3 + 21153A 4 ft 4 + ... = 0.(46) 

It is important to notice here that the non-perturbative character of the 
approximations studied in this work is due to non-perturbative treatment 
of the gap equation - working with an effective action at a given order in 
ft the gap equation will be truncated to that order in ft and solved after 
setting ft — 1. One can easily check that if one expanded the solution of 
the gap equation O to the given order in ft, the results for the ground state 
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energy and for the second excitation would coincide to that order in A with 
energies obtained in perturbation theory for Schrodinger equation. A field- 
theoretical derivation of the perturbative result for the second excitation in 
the loop expansion of the conventional effective action is considerably more 
laborious, since more diagrams have to be evaluated 0. 

In this approach approximations to higher excitations cannot be obtained, 
because composite propagators for higher number of particles cannot be de- 
rived directly from T[A]. There is however another strategy, analogous to 
Brillouin-Wigner perturbation theory, which seems to provide a way out, by 
calculating all zeros in the expression for T 4 (p) of this order, and treating 
them as approximations to successive excitations. We have checked that the 
lowest zero of T A (p) provides a good approximation to the second excitation, 
and in this case the result is close to (f45|) . However, the other zeros do not 
agree with the exact results: the higher the excitation, the worse the agree- 
ment. Even the perturbative behavior of the excitations, obtained by using 
a solution of the gap equation expanded to the given order in h, is wrong. In 
fact, the Brillouin-Wigner strategy applied to the two-particle inverse propa- 
gator is not able to provide reasonable approximations for higher excitations 
and does not show any supremacy over the Rayleigh-Schrodinger strategy. 
Therefore, in our work we use the Rayleigh-Schrodinger strategy, approxi- 
mating excitation energies by zeros of the composite propagators calculated 
to the considered order in h. In this approach higher excitations can be inver- 
stigated by considering the effective actions for higher composite operators 
which directly generate higher composite propagators. 



3.1.2 The approximations derived from T[A,A] 

We shall derive an approximation to the AO spectrum from the effective 
action for $ 2 (x) and $ 4 (x), T[A, A], discussed in Section ([2.3|). Calculating 



the Feynman diagrams of Fig. ?? in the space-time of one dimension the 
ground state energy is obtained in the form 

Eo = h^-(n 2 -m 2 )h^+3\h 2 A 2 +~A 2 h 3 n 5 +h 3 \A-4A 3 h*n 7 +..., (47) 
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where ( f4"3|) can be used to eliminate A in favor of Q. The inverse two-particle 
propagator, taken from Fig. ??, in one-dimensional space-time becomes 

r 4 ( p ) = -(p 2 + AQ 2 ) + 6\n-i6A 2 Q 7 h 2 160Q2 + p2 



2^ ' 3(160 2 + p 2 ) 

3 3 9 256n 4 + i76ny + p 4 
+ 64 ^ Afi (i 6 n* + p 2)2 — + (48) 

and its zero, calculated to order ft 3 , determines the second excitation to be 

„ 3Aft in AW 9A 2 ft 2 _ A 3,3^7 nnA 2 ,3^ 2a 27A 3 ft 3 
e 2 = 20 + _ - 104— _ - 96A W + 20AWA + _ + . 

The inverse four-particle propagator, from Fig. ??, takes the form 
r 8 G») = 20 3 (160 2 +;p 2 ) + ft(320 5 - 576A0 7 + 20y-12AOy) 

h 2 

+ — TTTo ow ^ — f414720 9 - 746496A0 11 - 331776A 2 13 

9(40 2 + p 2 )(360 2 +p 2 ) V 

+ 141120V - 222912AOV + 544384A 2 O n j9 2 + 1008Oy - 9504AO 
+ 32768A 2 OV + 18flV - 108AO 5 / + 88A 2 V 

and its zero, calculated to order ft 2 , determines the fourth excitation to be 

£4 = 4n-24A ft tf-^^! + ... (51) 

In the above expressions the effective mass, O, and the effective coupling, A, 
are space-time independent and have to be determined from the gap equa- 
tions (||) and (fjl), which in one- dimensional space-time become 

2 02 o*ao4 176ft 2 A 2 6 400ft 3 A 3 8 

m -0 2 -8ftA0 4 H h... = (52) 

3 9 

5fiS97) 3 A 3 9 

24ftA + 32ftA0 5 - 288ft 2 A 2 7 + — — ... = (53) 

When considering the approximations to the AO spectrum at a given order, 
we shall truncate these equations to that order in ft, and solve them after 
setting ft = 1. If one expanded instead the solutions of the gap equations, O 
and A, in powers of ft, the perturbative results for the ground state energy 
and the second excitation to order A 3 and for the fourth excitation to order 
A 2 would be recovered. 
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3.2 The numerical results for the AO energy levels 

Here we shall compare the numerical results for the AO spectrum, obtained 
from T[A], with those from T[A, A]. With the use of T[A] only the ground 
state and the second excitation can be calculated, but T[A, A] enables us to 
obtain in addition the fourth excitation. Both methods are non-perturbative, 
since the gap equations truncated to the given order in h have been solved 
numerically, after setting h = 1. In both cases the solution with the largest 
positive value for Q has been chosen. The results are compared with the 
exact values of the AO energy levels, calculated by a numerical procedure 
based on a modification of the linear variational method jl7|, and with the 



results of perturbation theory. All results are presented as functions of the 

2 

dimensionless quantity z = 7^373, which is the only parameter of the theory 
after rescaling all quantities in terms of A. 

In Fig. 10 we present the results for the ground state energy to successive 
orders in h. The quality of an approximation can be expressed in terms of 
a critical value of the parameter z, below which large discrepancies between 
results of different orders and the exact value appear. The critical value for 
approximations ([42]) obtained from T[A], z crit « — 1., appears much smaller 
than that for the perturbative results, z crit « 2.. The approximations (|47|) , 
obtained from T[A, A], are still better hence the critical value is still smaller, 

Zcrit ~ 3. 

The results for the second excitation, shown in Fig. 11, are very similar to 
those for the ground state, only the critical values are correspondingly larger: 
z cr it ~ 4. for the perturbative calculations, z cr %t ~ 2. for the approximations 
obtained from T[A] (|45"D, and z cr u ~ —.5 for that obtained from T[A, A] (??). 
These results should be compared with those obtained in the same order of 
the 2PPI expansion ||10|| . We observe that the two methods of expanding 
T[A] provide the approximations to the AO spectrum of a similar quality. 

The approximations for the fourth excitation can be obtained from the 
inverse four-particle propagator T 8 (p), which can be derived directly from 
T[A, A], but not from T[A]. As can be seen in Fig. 12 the results (|5T|), obtained 
from T[A,A], are in good agreement with the exact result for the values of 
the parameter z greater than z cr u ~ 1. This should be compared with the 
perturbative results where z crit ^5. 

It is interesting to observe that T[A] and T[A,A] provide reasonable re- 
sults for ground state energy even in the case of the double-well anharmonic 
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potential (z < 0), if the wells are not too deep. Since the results derived 
from T[A, A] are considerably better, one may suppose that including higher 
composite operator could improve the convergence properties and make the 
method applicable even in the case of deeper wells. In all methods consid- 
ered the numerical results for the excitations are worse than that for the 
ground state: the higher the excitation, the larger the critical values below 
which the method does not converge, but in any case including the operator 
$ 4 improves the convergence. One can speculate that increasing the num- 
ber of composite operators included in the effective action, we could make 
the method applicable for higher excitations in the case of the double well 
potential. 

4 Conclusions 

In the quantum scalar field theory the effective action for the local com- 
posite operator $ 2 (x), T[A], can be obtained by performing the Legendre 
transform of the connected generating functional W[J] order by order in h. 
This provides a convenient tool to study the ground state and the second 
excitation in the theory. The result to the order h 2 coincides with the Gaus- 
sian approximation, where the effective propagator is a Hartree one; higher 
orders provide a way to go beyond this result. A further improvement of 
the method can be achieved by taking into consideration higher composite 
operators <3> fc (x); their effects appear successively in higher orders of h. 

In this work we investigated the effects of including the operator $ 4 (x) 
in the effective action, T[A,A], which has been calculated as a series in h. 
The effects of the operator $ 4 (x) first appear at order h 3 . The successive 
approximations for the energy density and for the inverses of the two- and 
four-particle propagators, T A (p) and T 8 (p) have been derived from T[A,A]. 
To each order in h the effective mass and coupling have to be determined from 
the algebraic gap equations. The numerical results for the ground state and 
lowest excitations have been calculated for the theory in the space-time of one 
dimension, i.e., for quantum mechanical anharmonic oscillator. In this case 
the gap equations were solved numerically, and the excitation energies were 
determined as zeros of the inverse composite propagators in successive orders 
of h. A comparison with the results obtained from the effective action for 
the operator <3> 2 (x), T[A], shows that including the operator $ 4 (x) improves 
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the convergence of the approximation to the ground and the second excition 
considerably. Moreover, the effective action T[A, A] enables us to obtain the 
approximation to the fourth excitation which is not possible with the use of 

r[A]. 

One may note that the Gaussian approximation, obtained as the lowest 
approximation to the effective action for the operator $ 2 (x), appears also 
in a variational calculation with the Gaussian trial propagator. However, it 
is not staightforward to obtain a substantial improvement beyond the result 



of this simple ansatz by modification of the trial propagator fllSfl . Including 
higher operators into the effective action provides such a improvement, and 
in addition enables one to study excitations of the system. 
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Figure captions 



Figure 1. The loop expansion of W[J,K] in terms of Feynman diagrams in 
configuration space; the line denotes the propagator (— d 2 + Q 2 (x))~ 1 5(x — y), 
and the full circle stands for the four-particle coupling (— 24A + K(x))S(x — 
y)S(x — w)5(x — z). 

Figure 2. A [J, K] and A [J, K] represented with Feynman rules of Fig.l; the 
small empty circle denotes an external point x. 

Figure 3. J [A] obtained by inversion of A [J] in powers of h. The line de- 
notes the propagator (— d 2 + fl 2 (x))~ 1 5(x — y) and the four-particle vertex 
is — 24X5(x — y)5(x — w)5(x — z). Two lines meeting at two points with the 
slash across represent the inverse of the free two-particle propagator. 

Figure 4. T[A] represented with Feynman rules of Fig. 3. 

Figure 5. The inverse of the two-particle propagator, calculated from T[A], 
represented in terms of Feynman diagrams in momentum space; the line 
stands for a free propagator 2 + n 2 ; the small empty circle denotes an exter- 
nal momentum p. 

Figure 6. J[A, A] and K[A, A] obtained by inversion of A [J, K] and A[J, K] in 
powers of k; the single line denotes the propagator (— d 2 + m 2 — J (x))^ 1 5(x — 
y), the lines meeting at two points represent the free composite propagator 
of the corresponding number of fields, and the slash means the inversion. 

Figure 7. T[A, A] represented with Feynman rules of Fig.6. 

Figure 8. The inverse of the two-particle propagator, calculated from T[A, A], 
represented in terms of Feynman diagrams in momentum space; the line 
stands for a free propagator p2 ^2 , the small circle denotes an external mo- 
mentum p. 

Figure 9. As Fig. 8, but for the inverse of the four-particle propagator, cal- 
culated from r[A, A]. 
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Figure 10. The ground state energy of the AO, obtained to the given order 
of h from T[A] (dotted lines) and from T[A, A] (dotted- dashed lines), plotted 
vs z — ^573-; compared with the exact value (solid lines) and given order 
perturbative results (dashed lines). 

Figure 11. As in Fig. 10, but for the second excitation energy of the AO. 
Figure 12. As in Fig. 10, but for the fourth excitation energy of the AO. 
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Figure 1: The loop expansion of W[J,K] in terms of Feynman diagrams in 
configuration space; the line denotes the propagator (— <9 2 + m 2 — J(x))~ 1 6(x — 
y), and the full circle stands for the four-particle coupling (—24:X-\-K(x))6(x — 
y)S(x — w)6(x — z). 




Figure 2: A [J, K] and A [J, K] represented with Feynman rules of Fig.l; 
small empty circle denotes an external point x. 
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Figure 3: J[A] obtained by inversion of A [J] in powers of fi. The line de- 
notes the propagator (— d 2 + ft 2 (x))~ 1 S(x — y) and the four-particle vertex 
is — 24A<5(x — y)S(x — w)6(x — z). Two lines meating at two points with a 
slash across represent the inverse of the free two-particle propagator. 
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r[A] = J[|(0 2 (x) - m 2 )A(x) - 3Xh 2 A 2 (x)]dx 
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Figure 4: T[A] represented with Feynman rules of Fig. 
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Figure 5: The inverse of the two-particle propagator, calculated from T[A], 
represented in terms of Feynman diagrams in momentum space; the line 
stands for a free propagator 2 _^2 ■> the small empty circle denotes an external 
momentum p. 



J[A, A] = m 2 - n\x) + 12h\A(x) - §AK[A, A] 




Figure 6: J[A, A] and K[A, A] obtained by inversion of A [J, K] and A [J, K] 
in powers of %\ the single line denotes the propagator (— d 2 -\-Q 2 (x))~ 1 S(x — y) } 
the lines meeting at two points represent the free composite propagator of 
the corresponding number of fields, and the slash means the inversion. 
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r[A, A] = J[|(0 2 (x) - m T )A(x) - 3\h 2 A 2 (x) - h 3 \A(x)]dx 




Figure 7: T[A, A] represented with Feynman rules of Fig. 6. 
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Figure 8: The inverse of the two-particle propagator, calculated from T[A, A], 
represented in terms of Feynman diagrams in momentum space; the line 
stands for a free propagator p2 + n 2 , the small circle denotes an external mo- 
mentum p. 
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Figure 9: As Fig. 8, but for the inverse of the four-particle propagator, cal- 
culated from r[A, A]. 
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Figure 10 




Figure 11 
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